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We present a systematic theoretical analysis of the motion of a pair of straight counter-rotating 
vortex lines within a trapped Bose-Einstein condensate. We introduce the dynamical equations of 
motion, identify the associated conserved quantities, and illustrate the integrability of the ensuing 
dynamics. The system possesses a stationary equilibrium as a special case in a class of exact solutions 
that consist of rotating guiding-center equilibria about which the vortex lines execute periodic 
motion; thus, the generic two-vortex motion can be classified as quasi-periodic. We conclude with 
an analysis of the linear and nonlinear stability of these stationary and rotating equilibria. 



I. INTRODUCTION 

Vortices are persistent circulating flow patterns that 
occur in many diverse scientific and mathematical con- 
texts [1], ranging from hydrodynamics, superfiuids, and 
nonlinear optics to specific instantiations in sunspots, 
dust devils [2], and plant propulsion [3]. In atomic 
Bose-Einstein condensates (BECs) at ultracold tempera- 
tures [4-6], quantized vortices arise as especially persis- 
tent topological defects that play important roles in both 
Hamiltonian and dissipative dynamics, as well as in quan- 
tum turbulence. Vortices and vortex lattices in atomic 
BECs have been analyzed in a series of reviews [7-9]. 

Recently, much theoretical [10-17] and experimen- 
tal [18-21] attention has been devoted to small vortex 
clusters, with a special focus on the simplest "vortex 
molecule," the vortex dipolc. This structure, consisting 
of two countcrcirculating vortices, is central to the re- 
laxation of superfiuids through pairwisc rcconncction of 
the vortex lines. Such reconnections have recently been 
observed in turbulent superfluid helium [22], a system 
considerably more complicated than the dilute-gas BEC. 

Vortices in a Bose-Einstein condensate are observed 
experimentally by identification of the vortex lines (or 
"cores"), which are regions of reduced atomic density 
that surround the phase singularity about which the fluid 
rotates. The dynamics of the vortex lines in a dipolc 
are, in general, three-dimensional (3D), and can be sim- 
ulated through a mean-field model, namely the Gross- 
Pitaevskii partial differential equation (PDE) [18, 23- 
26]. In many cases of experimental interest the vortex 
lines arc straight and parallel, enabling calculations based 
on two-dimensional (2D) reductions of the PDEs. Such 
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simulations have not proven entirely satisfactory in their 
agreement with recent experimental results [19, 21]. An 
alternative 2D formulation is to focus on the dynamics 
of the vortex lines themselves by evaluating an ordinary 
differential equation (ODE) that treats the locations of 
the lines as individual (quasi) particles [16, 21]. The lat- 
ter approach privileges the vortex lines themselves and 
enables considerable insight into their fundamental dy- 
namics, even if the loss of the third dimension precludes 
a direct investigation of reconnection dynamics. 

Our aim in the present work is to offer a systematic 
study of the equations of motion ensuing from the quasi- 
particle approach, focusing on an ODE model of two 
straight counter-circulating vortex lines in a cylindrically 
symmetric Bose-Einstein condensate. We calculate the 
conserved quantities within the model and demonstrate 
its integrability in the case of two vortices. The vor- 
tex lines therefore undergo generically quasi-periodic mo- 
tion. To cement this observation, we present a system- 
atic analysis of the stationary or time-periodic states that 
the system possesses. We identify time-periodic guiding 
center solutions and explicitly compute their frequency 
as well as the fixed distances of their constituent vor- 
tex lines from the center of the superfluid. Finally, we 
examine both the linear and nonlinear stability of such 
states, arguing that the generic quasi-periodic motion of 
the two- vortex system consists of such guiding-center ro- 
tations along with (second frequency) epicycles around 
these states. 

Our presentation is structured as follows. In Section II, 
we present the mathematical model, while in Section III, 
we examine its conserved quantities, integrability and as- 
sociated dynamics. In Section IV we consider the special 
stationary equilibrium, which we generalize to rotating, 
guiding-center equilibria in Section V. Finally, we con- 
clude our presentation in Section VI with some interest- 
ing aspects of the problem for future study. 
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FIG. 1: Elements of the vortex particle model, (a) In an 
unbounded homogenous fluid, two counter-circulating vortex 
lines move together (dashed lines) in the direction of the flow 
between them, (b) A single vortex line in a bounded conden- 
sate precesses (dashed line) about the center of the condensate 
O in the same sense as its circulation. The radius R may be 
either the Thomas-Fermi radius of the condensate or the wall 
of the confining cylinder, (c) Two counter-circulating vortex 
lines in a bounded condensate experience instantaneous veloc- 
ities v''^^'^'' and v^"2 ' associated with precession and interac- 
tions, respectively (thick arrows). The thin arrows define the 
displacements of the two vortices with respect to one another 
(r2i) and the center of the condensate (ri,2). 



scopic precession about the condensate center at fixed 
distance rj (Fig. lb). The precession frequency may be 
expressed in both cases by [7, 9, 19] 
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where /i is the chemical potential, tOr and tOz are the 
confining radial and axial frequencies of the harmonic 
trap, respectively, and R is the Thomas-Fermi radius of 
the condensate or the radius of the hard- wall cylinder [9] . 
The common dependence on is what draws together 
these two different physical situations. 

Summing these two velocities for the kth vortex in a 
"gas" of n vortices gives, in complex coordinates [21], 



II. MODEL 

We consider straight line vortices in a Bose-Einstein 
condensate with cylindrical symmetry. The most promi- 
nent physical example arises in the context of harmon- 
ically confined, oblate BECs [27], but the same general 
equations and conclusions may be drawn in the case of a 
hard- wall container as well [9]. We begin by considering 
an ensemble of singly-quantized vortices, and in the fol- 
lowing section specialize to the case of the vortex dipole. 

The motion of a vortex line involves its interactions 
with other vortex lines, as well as the effect of "boundary 
conditions" associated with the fiuid confinement [16]. 
We begin with the interactions. The angular velocity of 
the fluid flow a distance r from a singly-quantized vortex 
core in a homogenous fiuid is 



-Sk^{rk)zk + ^3^i''3k){zk - Zj), (4) 
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where m is the atomic mass. This flow pattern deter- 
mines the motion of a second vortex line, which moves 
with the flow at its center (core) [28]. The motion of 
the flrst vortex line is similarly determined by the flow 
pattern established by the second. If the separation 
between the jth and kth. vortex lines is rjk, then the 
two same-charge vortex lines orbit one another at an- 
gular frequency ^{rjk)- On the other hand, if the vor- 
tices are counter-circulating (vortices of opposite charge) , 
then the vortex lines move together with linear speed 
Vjk = Tjk^irjk) in the direction of the flow between them 
(Fig. la). 

The motion of the kth vortex line in either a hard- wall 
or a disk-shaped harmonic potential also involves gyro- 



where {xk,yk) is the position of the kth vortex and 
Zk = Xk + ivk = ■rfee*^% Vk = \zk\, and r-jk = \zk - Zj\. 
The topological charge of the fcth vortex is Sk — ±1, 
with the positive (negative) sign referring to counter- 
clockwise (clockwise) circulation as viewed from the pos- 
itive z axis. The constant parameter 6 modifles the in- 
teraction strength slightly from the homogeneous case in 
a harmonic trap; for experimentally relevant parameter 
values, it has been argued that to a good approxima- 
tion it is given by 6 = 1.35 [16], while for the hard-wall 
potential it is 6 = 2. 

The relevant velocities and coordinates for a vortex 
dipole are shown in Fig. Ic. Since the model is only 
valid for straight line vortices, we adopt the convention 
of referring to the location of a vortex line as the location 
of the vortex (particle). 



III. CONSERVED QUANTITIES AND 
INTEGRABILITY 



The system of interest has a Hamiltonian structure. To 
see this, let us write the differential equations in Carte- 
sian coordinates as 

±k = -Sk^{rk)yk - 9 X] ^o^i'^jk)iyk - Vj), 



(5) 
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For simplicity we work henceforth in dimcnsionlcss dis- 
tance units with R = 1. Let us define the function 

O " 

ff(zi,...,z„) = -tfj2Hl-rl) 

fe=i 



cl + yl and = b^{R) = hb/m 
with R = 1. Then, it can readily be found that 



where Vk = \zk\ = ^/xt 
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for every k = which means that iJ is a first 

conserved quantity, i.e., the first integral or Hamiltonian 
along the orbits of the system. Given initial positions 
Zk{0), then 



Hiziit),...,z„{t))^Eo yt, 
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for some suitable constant Eq. 

It is directly verifiable that a second conserved quantity 

is 



krl, 
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which represents the angular momentum of the system; 
see, e.g., Ref. [29] for a discussion about conservation 
laws (in the absence of the precessional terms). 

The existence of two conserved quantities guarantees 
integrability in the classical Liouville sense [30] for the 
case n = 2, whether the vortices are co- or counter- 
rotating. This implies that the energy level sets are com- 
pact and the phase space is foliated by invariant tori. On 
each of these, the motion is quasi-periodic with two fre- 
quencies. In the following, we consider other dynamical 
aspects of the n = 2 case for vortex dipoles with = 1 
and 5*2 = —1. 



A. No Collisions 



that is, the vortices are separated by a computable min- 
imal distance that depends on the initial position of the 
vortices. 

It is worth reiterating that the model considers the 
dynamics of straight vortex lines only, such as those oc- 
curring in oblate BECs, and does not therefore preclude 
the possibility of collisions and possible reconnection phe- 
nomena when the vortex lines are tilted or bent. 



B. No Ejections 

The vortices never reach the edge of the fluid, remain- 
ing confined instead within a computable inner circle. 
Assuming ri2(t) < 2 for all t (i.e., the distance between 
the two vortices cannot exceed twice the radius, R = 1, 
of the condensate), we have 



(l-rf)''«(l-r^)''" 
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In consequence, 



rUt) < 1 



Vt. (12) 
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where h = (c22-*«/2)i/f2o < i 

Again, it is worth noting that at finite temperature 
the dissipative presence of thermal atoms is expected to 
cause the vortices to leave the fluid. These effects are not 
considered in the present model. 



IV. THE STATIONARY EQUILIBRIUM AND ITS 
STABILITY 



We now prove the existence of a stationary equilibrium 
of the vortex pair dynamics and illustrate its stability. 



A. Existence 

Lemma 1 There is an equilibrium, unique up to rota- 
tions, given by 
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The two vortices never collide. To see this, we exploit 
the fact that the Hamiltonian is constant along orbits. 
Taking exponentials on the Hamiltonian, we have 



(l-r?)"«(l-r2)"o^*«/2^C2 >o 



(10) 



where = (1 - n (0)2)^^« (1 - r2(0)2)f2ori2(0)*«/^ Note 
that < < 2*«/2 because < r,{0)^ < 1 (i = 1,2) 
and < ri2(0) < 2. A first consequence of Eq. (10) is 



ri2(t) > 



Vt, 



(11) 
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Proof. Recall that ?i = 2 and Si = -52 = 1. The 

system (4) in complex notation is 



izi = -n{ri) zi + -^^{zi - Z2) 

^R 

iz2 = +^{r2) Z2 - -z-2~{^2 - 21) 



(14) 
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Here Zi = ri exp{i9i) and Z2 = ^2 exp(i6'2) are, in general, 
time-dependent, but we look for equilibria, i.e., constant 
solutions. By direct substitution of the given solution 
above, one sees that it is indeed a constant solution. Let 
us prove that this is unique. Naturally, ri,0i denote polar 
coordinates (or equivalently an amplitude-phase decom- 
position). We can find that the evolution of such polar 
coordinates reads: 
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^' 12 
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Let us call 6 = 61 — 62. From the last equation, sin6' = 0. 
Then, 6 ^ or 9 = it. It is straightforward to show that 
the first option must be discarded. In particular, ii 6 = 0, 
from the first equation 



B. Stability 

Theorem 1 The equilibrium (r]',r2,0o) stable (in the 
sense of Lyapunov). 

Proof. In the coordinates (ri, r2, 6*), the system of equa- 
tions (14) reads 



ri 



r2 



n{ri)+n{r2) 

'^Rr2 
^'12 



2_!1 _H 

r2 ri 



cos 6, 
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Of course, H defined by Eq. (6) is still a conserved quan- 
tity; in the new variables, it reads 



H{ri,r2,e) 



1 
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\Ti{r\ + r2 — 2rir2 cos ( 



= n{ri) 



Since ^(ri) > 0, 
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Note that basic trigonometrical considerations give = 
r\+ r2 — 2rir2 cos 6. 

The equilibrium (r", r2, 6'o) is a critical point of H . The 
Hessian matrix evaluated at (r]',r2,6'o) is 



2r|, 



f^(ri) >0, 



and then ?'i > r2. Arguing analogously with the second 
equation, one gets r2 > ri and thus we reach a contradic- 
tion. In conclusion, it must be the case that 9 = n. Now, 
we prove that ri = r2. On the contrary, let us assume 
that ri > r2 (the remaining case is analogous). Adding 
the two first equations 



0„2 
^'12 



r2_ 

ri 
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r2 



<0, 



which leads to a contradiction because 57 (r) is strictly 
increasing. 

In conclusion, 6 = ir, ri = r2 and wc can find the exact 
value 7'i = r2 = 40^+^) ^ t>y solving the first equation 
of the system (15). 

To emphasize the rotational invariance of our system, 
it is convenient to work with the variables (ri ,r2,9). The 
(unique) equilibrium then reads 



V 4r!o + $i?'' 
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One can easily prove that this matrix is positive-definite 
by Sylvester's criterion. Hence, iJ is a Lyapunov function 
that attains a minimum at (rj,r2,0o)j and the equilib- 
rium is stable. 

The stationary equilibrium is relevant to vortices that 
are symmetric with respect to a line passing through the 
condensate center. We illustrate several typical orbits in 
Fig. 2. This class of orbits has been predicted theoreti- 
cally [15, 16] and observed experimentally [18, 21]. For 
relatively small perturbations from the symmetric equi- 
librium, the (linear) rotational frequency (i.e., the imag- 
inary part of the stability eigenvalue associated with the 
equilibrium configuration) of these orbits is in reasonable 
agreement with those found experimentally [21]. 



V. GUIDING CENTER EQUILIBRIA AND 
THEIR STABILITY 

Finally, we prove the existence of a class of guiding 
center equilibria and illustrate their stability. 
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FIG. 2; (Color online) Typical computed orbits associated 
with symmetric vortex configurations, determined by inte- 
grating the equations of motion [Eqs. (14)]. The vortex on 
the left (blue) has a counterclockwise sense of circulation; the 
vortex on the right (red) has a clockwise sense of circula- 
tion. The stationary equilibrium positions of the vortices are 
indicated by stars. These orbits arise through competition 
between the effects of interaction and precession [15]. 



A. Existence 

Consider two vortex positions {xi,yi) and {x2,y2), 
or in complex notation, zi = ri exp(i0i) and Z2 = 
r2cxp(i6'2), where all of the r's and 9's are time- 
dependent. The system under study is 
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— 2-(^2 - Zi). 
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(17) 
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We seek precessing guiding centers about which the 
vortices oscillate. For the special case of guiding cen- 
ters equidistant from the center of the fluid, we expect 
the precession frequency uj to vanish, reproducing the 
stationary equilibrium case above. Given a vortex posi- 
tioned at one guiding center, and a (counter-circulating) 
anti-vortex at the other, we expect that each vortex will 
stay at its own guiding center; that is, if the guiding 
centers are themselves precessing, each vortex orbits the 
center of the condensate at the precession frequency of 
the guiding centers. From symmetry, we also expect the 
angular positions of the guiding centers to be tt radians 
apart on a line that passes through the center of the con- 
densate, i.e., on opposite sides of the center. 

Therefore, we adopt a trial solution of the form 

Zi{t) ~ riexp(ia;i), 

Z2{'t) — ^2 exp(ia;t + ivr) = — 7'2 exp(ici;t). 



where ri and r2 are now constant, as is w, the preces- 
sion frequency of the guiding center. Note that in this 
case ri2 = ri + r2 = s, the constant separation distance 
between the two guiding centers. 

It is now convenient to introduce the notation: 



1 



/3 = 



1 



and 



1 



7 = 



^' 12 
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with a, j3, and 7 all time- independent constants, by as- 
sumption. 

Since the time dependence in zi{t) and Z2{t) is now 
explicit, we take the derivatives directly: 

zi = zwri exp(iwt) and 22 = — iwr2 exp(iwt), 

whereupon the differential equations become 

wri = af^o^'i - 7*i?.(?'i + r2), (19) 



ur2 



(20) 



Recalling s = ri + r2 = (27)^^/^ as a final (as yet un- 
used) constraint, we now have three equations and three 
unknowns (ri, r2, and w); that is, upon choosing s, the 
guiding center locations ri and r2 can be identified in 
terms of s, as well as the precession frequency of the 
guiding centers uj — provided that there are solutions to 
this set of equations. Given the nature of the ensuing al- 
gebraic equations (see below), such solutions will generi- 
cally exist. We derive the relevant equations below. 

We start with the precession frequency u. From the 
system of Eqs. (19) and (20), we find 
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Note that if r2 = ri, then a ~ l3 and cj = 0, recovering 
the expected stationary equilibrium. 

Furthermore, from the same Eqs. (19) and (20) we im- 
mediately obtain: 



17o(a + /3)-7$i^ (2 + ^ + 

7-2 
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0. 



(22) 



Substituting in for a, /?, and 7, in terms of ri and r2, 
and simplifying, one finds 



no 
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Once the position of r2 is fixed, the location of the first 
vortex r = ri can be found from the third order polyno- 
mial 



I3r' - ~ (no + P)r + ^ = 0. 

2^2 7-2 



(24) 



An alternative possibility is to fix s, then, the location 
of the first vortex r = ri can be found from the quartic 
algebraic equation 



C4r^ + c^r^ + C2r^ + Cir -f- cq = 0, 



(25) 
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(26) 



where the coefRcients Cj (j = 0, 1, 2, 3) are given by: 

CO - B(s'-l), 

ci = 2s- - 2sB, 

C2 = 3.5^ - 2 - B{s^ - 2), 

C3 = -4.s + 2sS, 
C4 = 

and B = $fl/(217o). 



B. Stability 

In the following we prove the stability of the guiding 
center of the form 

zi(i) = ri cxp(iajt), Z2(t) = — ^2 exp(iajt), (27) 

obtained in the previous subsection. To this purpose, 
wc pass to a co-rotating frame by making the change of 
variables Zi = Zi cxp{~iujt) on the original system (14). 
Keeping the more convenient Zi notation for the state 
variables in a slight abuse of notation, the resulting sys- 
tem is 



IZi 



LUZi - Vt{ri) zi + -^^-^{zi - 22), 



12:2 = ^^Z2 +Q{r2) Z2 



^'12 
^' \2 



(28) 



{Z2 - 21). 



Note that a rotating solution like (27) of the original 
system (14) is equivalent to the statement that (y'l, — r2) 
is an equilibrium of the new system (28). 

In Cartesian coordinates, the system (28) reads 



^'12 



±2 = ujy2 + ^{r2) y2 - *_R. 
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~ol2 ' 
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y2 - yi 
^'12 
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and passing to polar coordinates, we obtain 
= _a;-|-f7(ri) 

62 = - Vt{r2) + 
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ri = „ 2 ^in(^i - ^2), 
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— cos(t^i — 
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(30) 



r2 



2r^2 
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2— sin(6'i - 6I2). 



Of course, when a; = we recover Eqs. (15). 
Now the Hamiltonian is 

H{ri,r2,e) = -'^{rl-rl) 



(31) 



ln(r^ + ^"2 ^ 2rir2 cos( 



where again 6 — 9i — 62- But here the additional diffi- 
culty is the absence of a simple explicit expression for the 
solution. 

Lemma 2 Any equilibrium of system (30) satisfies 9 — 

TT. 

Proof. From the last equation, 9 — -k oi 6 = Q. By 
contradiction, let us assume that 9 — 0. Then from the 
two first equations of system (30) 



2r2 
^' 12 



-w - n{r2) + 
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subtracting one from the other leads to 



n{n) + nir2) 
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2-H - H 
n r2 



0, 
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(32) 



then 



2-^-!:i>o 

ri r2 



and multiplying by rir2 one gets — (ri + r2)^ > 0, which 
is a contradiction. 

Now it is ensured that any equilibrium should be of 
the form (ri,r2,7r). We can then compute the Hessian 
matrix of the Hamiltonian above, evaluated at the point 
(ri , r2, tt) and obtain 

Q.{ri) + 5i + 5i2 - uj 5i2 

5\2 VL{r2) + 82 + 5x2 + 

(5i2rir2 



where 5, = -^^^ > and S12 - 2{ri+r2y 
is difficult to see that the matrix is positive-definite, but, 
from Eqs. (32), the frequency uj satisfies 



> 0. Here it 



oj = n{n) - <5i2 (^1 + ^ j = -n{r2) + 612 (^1 + ^ 

Inserting this information into the Hessian yields 
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FIG. 3: (Color online) Typical computed orbits associated with asymmetric vortex configurations, determined by integrating the 
equations of motion [Eqs. (14)]. (a) Motion of the vortices in the lab frame, where the outer vortex (blue) has a counterclockwise 
sense of circulation and the inner vortex (red) has a clockwise sense of circulation. The guiding center equilibria are indicated 
by stars, and their trajectories follow circles in the same sense as the rotation of the outermost vortex (black dashed lines), (b) 
Motion of vortices in a frame co-rotating with the guiding centers. Several representative orbits with the same guiding centers 
are shown, with the circle symbols attached to the curves corresponding to the motion in (a). 



and now one easily realizes that this matrix is positive- 
definite again by Sylvester's criterion. Hence, once again 
the relevant guiding-center equilibrium is a stable one, 
just as it is in the special case in which cj = and ri ~ r2. 

Rotating guiding center equilibria are relevant when 
the vortices are asymmetrically located with respect to a 
Hne through the condensate center. A representative ex- 
ample is given in Fig. 3; such behavior has been observed 
experimentally [21]. 

VI. CONCLUSIONS 

From the preceding calculations, we conclude that the 
vortex dipole system is an integrable one (this is true 
both in the co-rotating case, and in the counter-rotating 
one of principal interest herein). This implies that the 
energy level sets are therefore compact, and the phase 
space is foliated by invariant tori, on each of which the 
motion is quasi-periodic with two frequencies. In the 
context of the vortex dipoles, the first one of these fre- 
quencies is the precession of the guiding center equilibria 
which have been identified herein through the solution of 
relevant algebraic equations. The positive definite char- 
acter of the Hessian of the linearization in the appropriate 
variables (ri, r2 and 6 = di —62) guarantees stability of 
these equilibrium points. The second frequency is the 
oscillation frequency about the precessing equilibria. 

Our analysis and conclusions offer a straightforward 
view of the dynamics of vortex dipoles when the vortex 
lines are straight, such as those arising in oblate, harmon- 
ically trapped Bose-Einstein condensates. Nevertheless, 
many interesting questions arise. On one hand, in the 



counter-rotating vortex case, it seems particularly inter- 
esting and relevant to extend the considerations above to 
the case of three- or more vortex states within the con- 
densate and investigate the ensuing stationary [16] (and 
perhaps also guiding center) equilibria. On the other 
hand, another direction that naturally emerges concerns 
the examination of co-rotating vortices. In the latter 
case, we certainly expect the generalization of guiding 
center orbits, due to the common direction of rotation. 
Understanding the latter phenomenology in the general 
case of n vortices would arguably be interesting in its 
own right. 
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